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Abstract 

Equilibrium programming is investigated as a frame- 
work within which optimal structural design problems 
can be decomposed into subproblems. Aspects of equi- 
librium programming theory pertinent to this appli- 
cation are described, and three decompositions of a 
structural design problem into structural-sizing and 
structural-response subproblems are developed. Two 
of the decompositions are shown to be equivalent to ex- 
isting solution methods. However, the third method is 
novel and is shown to give optimal solutions both the- 
oretically, and in two numerical examples. 

Introduction 

Nonlinear mathematical programming (NLP) is ex- 
tensively used for optimal structural sizing in structural 
design. Utilizing NLP, a sizing design problem mini- 
mizes an objective function, such as weight, and satis- 
fies a set of design constraints, such as minimum gauge 
and local stress limits, by varying a set of structural- 
sizing design variables, such as facesheet thicknesses, 
and stiffener heights. The structural analysis that is 
required to compute the design constraints is subordi- 
nate to the optimization algorithm in the simplest use 
of NLP. That is, the structural response is recalculated 
for every change in the values of structural-sizing de- 
sign variables in the design process. The disadvantage 
of this approach is the expense of performing numer- 
ous structural- response and sensitivity- derivative calcu- 
lations during the design process. The incorporation of 
approximation concepts by Schmit and Miura (ref. 1) 
within NLP structural-sizing design methods signifi- 
cantly improved the computational efficiency of the 
method. Thus, approximation-based NLP structural- 
sizing design methods are now well developed, and are 
implemented in several commercial structural analysis 
codes. 

Though approximation-based NLP structural-sizing 
design methods have a computational advantage over 
simple NLP approaches, the design of structures incor- 
porating a high degree of detail requires a large number 
of design variables and constraints. Additional improve- 
ments in optimal structural- si zing methodology are re- 
quired to deal with increasingly large numbers of design 
variables and constraints. One approach to the solution 
of large problems is to decompose a large problem into 
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a set of smaller subproblems. Most of the smaller sub- 
problems can be solved separately and in parallel to re- 
duce the time required to find the optimum design. By 
repeating this decomposition process at the subproblem 
level, a multilevel decomposition (ref. 2 and references 
therein) can be obtained. One concern in the decom- 
position of an optimum design problem is whether any 
solution exists to the decomposed problem. If a solution 
exists, another concern is whether this solution to the 
decomposed problem is also a solution to the original 
problem. Although progress has been made in address- 
ing these concerns (see ref. 3), and in utilizing multilevel 
decomposition for optimal design, a mathematical the- 
ory that provides a framework for performing design 
decomposition appears lacking. 

The theory of equilibrium programming (EP) pro- 
vides a framework to analyze multiple, interdependent 
nonlinear programming problems. Equilibrium pro- 
gramming was developed in an operations research set- 
ting, and has been applied to economics, game theory, 
and network theory (ref. 4). Approaches to decompo- 
sition of structural-sizing design problems utilizing an 
equilibrium programming framework are developed in 
the present study. The decompositions examined in- 
vestigate several ways to separate the design problem 
into a structural-sizing subproblem, and one or more 
structural- response subproblems. The theorems on the 
existence of, and the necessary conditions for an EP so- 
lution are utilized to develop decompositions that have 
solutions, and to determine whether a solution is op- 
timal. In subsequent sections of the present paper, 
the characteristics of EP solutions are described, and 
some EP-based approaches to decomposition in struc- 
tural design are summarized. Specifically, structural- 
response and structural-sizing problems for minimum 
weight design in the presence of structural-sizing de- 
sign variable constraints, local buckling, stress, and dis- 
placement constraints are presented, and decomposed 
into several formulations utilizing EP. It is shown herein 
that two established structural- design methods; namely, 
fully stressed design (ref. 5) and NLP-based structural 
sizing using first-order Taylor series approximations of 
the constraints (ref. 1), can be derived using EP. A new 
decomposition derived using the EP framework is also 
presented. Computational results from this new EP- 
based formulation are compared with results from an 
established design method. 


Nomenclature 

Lower case symbols typed in boldface denote col- 
umn vectors, however some vectors are not denoted 
in this manner. The notation df /da for an arbitrary 
scalar / and a vector a = (ai,...,On) r denotes the 
gradient of scalar / with respect to vector a. This 
gradient is expressed as a row vector (i.e., df /da = 
(df/da\ t . . . ,dffdan))> The notation df/da for arbi- 
trary vectors f = (/j, . . . , fm) T and a = (a*, . . . , On) T 
denotes the gradient of vector f with respect to vector 
a. This gradient is expressed as an m x n matrix in 
which the component in column j of row t is given by 
dfi/daj . The notation df/da is utilized to represent nu- 
merical values of the total derivative of f with respect 
to the vector a using a matrix format similar to that of 
df/da . 

A vector of cross-sectional areas 

g vector of inequality constraint functions 

vector of external nodal forces 

/ objective function in nonlinear 

programming 

h vector of equality constraint functions 

K global stiffness matrix for a structural- 

response problem 

L Lagrangian function, (see equation (20)) 

M number of structural- response cases 

N vector of stress resultants within every 

element of the structure 

u vector of nodal displacements from a 

structural- response problem 

v vector of design variables for nonlinear 

programming, or vector of structural- 
sizing design variables in equilib- 
rium programming structural-sizing 
subproblems 

V strain energy of a structure 

W weight of a structure 

x, x 1 , x l vectors of all design variables for equi- 
librium programming, design variables 
of subproblem t, and all design vari- 
ables except those from subproblem i, 
respectively 

A vector of Lagrange multipliers for 

inequality constraints 

/i vector of Lagrange multipliers for 

equality constraints, (used for inequality 
constraints in the appendix) 


<y vector of stresses within every element of 

a structure 

Subscripts and Subscripts: 

A denotes an approximation of a function 

i denotes Pth equilibrium programming 

subproblem, or Pth load condition for 
structural response 

max maximum 

min minimum 

N denotes functional dependence on stress 

resultants 

s side constraints 

T transpose of a vector or matrix 

u denotes displacement constraints, or 

functional dependence on displacements 

cr denotes stress and local buckling con- 

straints, or functional dependence on 
stresses 

* denotes an equilibrium programming 

solution (i.e., the equilibrium point), an 
optimal solution, active constraints, or 
Lagrange multipliers corresponding to 
active constraints 

Theory 

To have a basis for comparison with the equilibrium 
programming (EP)-based formulations that follow, this 
section begins with a summary of the methodology 
for utilizing nonlinear programming (NLP) for struc- 
tural design. Some results from EP theory are then 
presented, and some characteristics of equilibrium pro- 
gramming problems are described. Three decomposi- 
tions of structural-design problems are developed using 
the EP framework, and some implementation issues are 
addressed. 

Nonlinear Programming Approach to 
Structural Design 

Nonlinear programming can be personified as having 
one decision maker (which may be implemented as 
a search algorithm) with control of a set of design 
variables given by the vector v. The goal of the 
decision maker is to minimize an objective function /(v) 
while satisfying a set of constraints. The mathematical 
description of a NLP problem is: 

min /(v) 

V 

g(v) > 0 0) 

h(v) = 0 
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Figure 1. Block diagram illustrating the major components 
of a solution algorithm for a NLP problem. 


where g(v) are inequality constraints - which could in- 
clude simple bounds on the design variables - and h(v) 
are equality constraints. A simple block diagram illus- 
trating the major components of a solution algorithm 
for a NLP problem is shown in figure 1. As shown in 
this figure, the analysis module is subordinate to the al- 
gorithm that searches the design space for the optimum 
design variables, and it can be queried for function (and 
gradient) values at will. 

The first-order conditions that are necessary for v = 
v* to be a solution to the NLP problem represented by 
statement (1) are often used in solution procedures. If 
a constraint qualification - such as linear independence 
of the active constraints at v = v* - is satisfied, then 
a set of Lagrange multipliers (A, p) exists such that at 
v = v* the following necessary condition relations are 
satisfied: 


d/(v) _ xT g g( v ) _ T gh ( y ) 
dv dv M dv 


= 0 


g(v) > 0 
h(v) = 0 


A > 0 


A T g(v) = 0 


( 2 ) 


The functions g(v) and h(v) for a structural- design 
problem implicitly involve the structural response; so 
results from structural analyses are necessary to deter- 
mine these constraint functions. 

Structural Analysis. A minimum potential energy 
formulation can be used to calculate the structural 
displacements and internal stresses that are required to 
evaluate a design. Given the structural arrangement, 


sizes for all the structural elements, a discretization of 
the structure into finite elements, and a set of external 
forces on the discretized structure F^, the correct 
structural displacements are those that minimize the 
potential energy of the structure. Thus, the structural 
response is the solution to the unconstrained NLP 
problem given by 


min (v^(u) - f £ u ) (3) 

where u is a vector of nodal displacements, V{n) is 
the strain energy of the structure, and F E is a vector 
of external nodal forces. In the present study, a lin- 
ear structural analysis is assumed, the strain energy is 
given by V(u) = l/2u r Ku, and the necessary condi- 
tion relations (2) for the unconstrained minimization 
problem represented by statement (3) simplify to the 
linear equations: 

Ku = F e (4) 

where K is the global stiffness matrix. Once the 
displacements u have been determined from equation 
(4), stresses <r within every element of the structure 
can be calculated from the displacements, the element 
strain- displacement relations, and the element elastic 
constants. 

Optimal Structural Sizing. A common objective 
function for structural sizing is the weight of the struc- 
ture W(v). The structural-sizing design variables v, 
referred to herein as sizing variables, can be the di- 
mensions of the individual elements that explicitly con- 
tribute to weight, such as beam dimensions, facesheet 
thicknesses, stiffener dimensions and spacing, and the 
thicknesses of a non-load-bearing spacer material; or 
they can be factors which effect the weight in an 
indirect manner, such as the orientation of compos- 
ite fibers. The constraint functions considered in the 
present study are the sizing variable constraints, or side 
constraints, g*(v); the local buckling and stress con- 
straints g£.(v,<r*(v)); and the displacement constraints 
gJi(u l (v)). The superscript i indicates that these quan- 
tities correspond to structural load condition i (i.e., to 
the vector of external nodal forces F^). The displace- 
ment constraint functions gu(u’(v)) are assumed to 
have no direct dependence on the sizing variables v, but 
the functional form u*(v) indicates an indirect depen- 
dence on v through the structural analysis. The stresses 
o'(v), which are shown to depend on v in the constraint 
functions, can have several forms. For example, one 
can write <r(v) = cr u (v,u(v)) = 0 >j(v,N(v, u(v))) to 
show that the functional form for stresses can depend 
directly on displacements u, or indirectly on displace- 
ments through the stress resultants N. Very often the 
side constraints, the displacement constraints, and the 
stress constraints reduce to simple bounds on the sizing 
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variables, the displacements, and the stresses, respec- 
tively. Thus, the NLP-based design approach in which 
the structural analyses are treated as subordinate to the 
sizing variable search algorithm can be summarized by: 


min W(v) 
g *(v) > o 
gu(u’(v)) > 0 
g^(v,<r’(v)) > 0 


( 5 ) 


where i = 1, . . . , M for M load conditions. The neces- 
sary conditions for v = v* to be the solution to the NLP 
structural- design problem given above are as follows: 


<W(v) ,T 9 g*( v ) 

dv 9 dv 


M 






chi* 


d\ 


+ (Kr 




dg £(v, er k ) dg£(v,q- fc ) d<r k {y) 

dv d<r k dv 

gs(v) > 0 


)] =° 


gu(u’(v)) > o 
gt(v,ff‘(v))>0 
A s >0 


A‘ u >0 

At>0 

(Aj) T gs(v) = 0 
(Au) T gu( ul ( v )) = 0 
(At) T gj r (v, < r i (v)) = 0 


( 6 ) 


where i = The terms u*(v) and du l (v)/dv 

are found by solving the following structural- analysis 
and sensitivity-derivative equations in which the param- 
eter vector v is assumed known: 


K(v)u* - F' E = 0 

k < v >£ + t-m 


(?) 


Sub- 

problem 



Figure 2. Block diagram illustrating the subproblem cou- 
pling in an EP problem. 


solution algorithms are summarized. References are 
made to the application of EP to design decomposition 
where appropriate. 

Equilibrium programming is a generalization of NLP 
which can be personified as having Af + 1 decision mak- 
ers (which may be implemented as search algorithms) 
that interact in a system. Each decision maker has a 
NLP subproblem to solve, and an independent set of 
design variables to control. To simplify notation, the 
design variables controlled by decision maker i are de- 
noted as x l , the design variables of all M + 1 decision 
makers are denoted as x = (x°, , . . , x^), and all design 
variables not controlled by decision maker i are denoted 
as x l = (x°,...,x l-1 ,x J+l ,...,x Af ). Decision maker i 
has an objective function to minimize, f i (x i ,5L i } } while 
satisfying a set of constraints. Thus, the mathematical 
description of equilibrium programming is: 

min /*(x*,x*) 
x* 

g’(x’,x*) > 0 ( 8 ) 

h*(x*,x’) = 0 


where it is also assumed that the nodal forces F^ are 
independent of v. 

Equilibrium Programming Theory 

In this section of the present paper, the formulation 
of an equilibrium programming problem is given. In 
addition, the first-order necessary conditions that are 
satisfied by a solution to an EP problem, a theorem on 
solution existence, some solution properties, and some 


for the t = 0, , . . , M interacting subproblems. The 
variables following the comma in any of the functions 
in statement (8) are treated as fixed parameters jn that 
subproblem. Thus, in the NLP problem of decision 
maker i, the design variables from other decision makers 
x* enter as parameters. The coupling of the M + 1 NLP 
subproblems through their design variables is illustrated 
in the block diagram of figure 2. In this figure, each 
subproblem is a NLP problem which might be solved 
by a solution algorithm like that of figure 1 . 
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The solution to all the NLP subproblems repre- 
sented by statement (8) is x = x*, which is called an 
equilibrium point. In a manner similar to NLP, first- 
order necessary conditions are satisfied at the equilib- 
rium point subject to a constraint qualification. Thus, 
at the equilibrium point x = x*, there exist Lagrange 
multipliers (A',/i l ) such that the following conditions 
are satisfied: 


a/V,**) ^T&i**,**) 

dx i [ } dx' 

iT ah^(x ^) =0 

^ ’ dx' 

gV,*‘)>0 ( 9 ) 

h’(x*,x’) = 0 
> 0 

(A*)V(x i ,x i ) = 0 


for i = 0, . . . , M . 

Existence and Optimality of an Equilibrium Point . 
In utilizing EP for a design decomposition, it is impor- 
tant that a solution exists for the decomposed problem. 
In fact, the requirements for existence of an equilib- 
rium point can be used as a guide in the decomposition 
process. A second concern involves the optimality of 
the equilibrium point. That is, is the equilibrium point 
resulting from the decomposition of a design problem 
also an optimal solution of the original NLP formula- 
tion? This second concern can be addressed by show- 
ing that EP necessary condition relations (9) for the 
decomposed problem imply the NLP necessary condi- 
tion relations (6) of the undecomposed problem. Both 
the existence and the optimality of an equilibrium point 
depend on the satisfaction of a constraint qualification. 

The satisfaction of a constraint qualification is re- 
quired for both the existence of a solution to the first- 
order conditions represented by statement (9), and the 
existence of an equilibrium point. One form for the 
constraint qualification is given in reference 6. It is sat- 
isfied if for all x feasible to the NLP subproblems repre- 
sented by statement (8), and for every i: 1) the vectors 
dhj(x % ,x l )/dx l are linearly independent, and 2) there 

is at least one solution z l to the relations 


dx 1 

3h i (x i ,x < )_ < _ n 

ft? Z 


( 10 ) 


where g* is the vector of inequality constraints in sub- 
problem i that are active at x. With regard to the in- 
equality constraints, this constraint qualification essen- 
tially states that it is always possible to move into the 


interior of the feasible region from a point on the bound- 
ary of that region. However, because constraint qualifi- 
cation relations (10) must be satisfied individually by 
the EP subproblems, the constraint functions in EP 
must satisfy more restrictive requirements than those 
in NLP. For example, the constraints from two EP sub- 
problems given by gj(-, x 2 ) and g 2 (x 2 , •) will not satisfy 
relation (10) individually because (dgl/dx l )z l s 0. 
However, relation (10) may be satisfied for this exam- 
ple when the constraints and design variables of the 
subproblems are combined within a single NLP prob- 
lem (i.e., ( dg,/dx)x > 0 where g. = (gj,g 2 ) and 
x = (x^x 2 )). 

A very general theorem for existence of an equilib- 
rium point, given in reference 6, requires continuity, 
but not differentiability, of the objective and constraint 
functions. Other conditions for existence are: 1) the 
functions satisfy constraint qualification relations (10) 
(actually only a weakened form of the constraint qual- 
ification is required), 2) the feasible region is bounded, 
3) the functions f l (x\x l ) and 0j(x\X*) are concave 
in x\ and 4) the functions /^(x^x*) are linear in x*. 
Other existence theorems are available which relax the 
concavity and linearity restrictions on (/\g*) and h x , 
respectively, if additional differentiability requirements 
are satisfied. 

Three methods are used in the present study for the 
decomposition of a design problem into EP structural- 
sizing and structural- response subproblems. The re- 
quirements for existence and optimality of the EP solu- 
tions are utilized as guidelines for applying these meth- 
ods in the decomposition process. In the first decom- 
position method, a constraint is applied only within a 
subproblem for which the design variables are effective 
in satisfying the constraint. Here effectiveness is de- 
termined by ability to satisfy constraint qualification 
relations (10). This method is akin to the concurrent 
subspace optimization method for optimization decom- 
position described in reference 7. In the second decom- 
position method, a constraint is applied within a cho- 
sen subproblem, and a change of variables is performed 
to yield a set of the design variables and constraints 
which satisfy the constraint qualification relations. In 
the third decomposition method, a constraint is applied 
within a chosen subproblem, and parameters of the con- 
straint functions which are normally independent of the 
subproblem design variables are replaced with approxi- 
mate models that depend explicitly on the subproblem 
design variables, thus satisfying constraint qualification 
relations (10). For example, in a constraint function 
g x (x x ) for subproblem t, the parameters x x which are 
normally independent of the subproblem design vari- 
ables x l are replaced with approximations that depend 
explicitly on x\ This method is closely related to the 
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approximation concepts of reference 1, and in fact, the 
use of explicit approximations of constraint functions in 
the NLP-based design methods can be viewed as a type 
of design decomposition. Examples of utilizing these de- 
composition methods in the development of structural- 
design equilibrium programming methods (SDEP) are 
given in a subsequent section of the present paper. 

Some Equilibrium Point Properties . Although the 
differences between an equilibrium point and an optimal 
point may appear slight, they are important. The fol- 
lowing equilibrium point properties, summarized from 
reference 4, illustrate the differences. Several examples 
of these differences can be found in the reference. 

An equilibrium point is, in general, different from 
an optimal point (i.e,, the solution of a NLP prob- 
lem), even if the same constraints are satisfied and each 
EP subproblem has the same objective function. This 
difference can occur because the coupling of the con- 
straint derivatives in the respective necessary condi- 
tions is generally less for an EP formulation than for 
a NLP formulation. Thus, the subproblems in EP can 
be called “loosely coupled.” Larger differences between 
the EP and NLP solutions can occur when the Af + 1 
subproblems in EP have Af + 1 different, and possi- 
bly conflicting, objective functions. (However, two of 
the structural- design decompositions described subse- 
quently are specifically formulated to ensure that an 
equilibrium point is also optimal.) An equilibrium point 
has a “stability” property which states that a solution 
to the EP problem does not change for a perturbation of 
the design variables of a single subproblem from equi- 
librium values. Additional constraints can affect EP 
solutions differently than NLP solutions. In NLP, ad- 
ditional constraints generally increase the value of the 
objective function. However, in EP it is possible for 
additional constraints to force a “cooperation” that re- 
duces the objective functions of all the EP subproblems. 

EP So/ut/on Algorithms . The solution of an EP 
problem can be obtained in several ways. The most 
straightforward method is to solve sequentially all 
the subproblems in some predetermined cyclical or- 
der. When the solutions to all the subproblems do not 
change from the previous iteration (within some numer- 
ical tolerance), the equilibrium point has been reached. 
This method is utilized in the present study. Other 
methods, such as solution of the set of nonlinear neces- 
sary condition equations (ref. 4), are also possible, and 
may be advantageous when convergence is difficult to 
achieve by the sequential solution method. 

Equilibrium Programming Approaches to 
Structural Design 

In this section of the present paper, decompositions 
of a structural- design problem into structural- response 


and structural-sizing subproblems are described. These 
decompositions are denoted as structural- design equi- 
librium programming (SDEP) formulations. In the suc- 
ceeding formulations, it is assumed that there are Af 
structural load conditions given by F l g } i = 1, . . . , Af ; 
and that, in general, the constraints of interest are siz- 
ing variable, stress, local buckling, and displacement 
constraints. In some cases, the formulations simplify 
appreciably when these constraints are reduced to sim- 
ple bounds. The formal definition of x 1 will be given 
for each subproblem. This formal definition may in- 
clude both optimization design variables, and behavior 
variables which result from a post- processing analysis 
within the subproblem. The conversion of the analy- 
sis for behavior variables to subproblem equality con- 
straints that conform to the EP formalism is obvious, 
and omitted in the development of the formulations. 

The three SDEP formulations that follow are de- 
rived using the methods for decomposition described 
previously. The logical steps for utilizing the decompo- 
sition methods to derive the formulations are described 
separately for each formulation. The first structural- 
design equilibrium programming formulation developed 
considers only stress, local buckling, and minimum 
gauge constraints. This method, termed SDEP1, is 
shown to be equivalent to the method of fully stressed 
design for rod and membrane elements. The second and 
third SDEP formulations, denoted SDEP2 and SDEP3 
respectively, provide for optimal designs with sizing 
variable, stress, local buckling, and displacement con- 
straints. Formulation SDEP2 is shown to be equivalent 
to the NLP-based approach to optimal structural siz- 
ing using a first-order Taylor approximation for a rapid 
analysis. Formulation SDEP3 is new, and the existence 
and optimality of its equilibrium points are discussed. 

Structured- Design Formu/ation SDEP1. It seems 
natural to define EP subproblems i = 1, . . . , Af , corre- 
sponding to the Af load conditions F^, to be structural- 
response subproblems, each a minimization problem 
represented by statement (3), and to define the subprob- 
lem design variables x* to be the displacements u l ob- 
tained from these minimization problems. This choice 
leaves subproblem 0 as the structural-sizing subproblem 
(i.e., minimize the structural weight by varying sizing 
variables x° = v subject to the sizing variable, stress, 
and buckling limits in the NLP problem represented by 
statement (5)). The shortcoming of this EP formula- 
tion is that an equilibrium point may not exist because 
constraint qualification relations (10) cannot be satis- 
fied for some constraints. For example, maximum stress 
constraints for rod elements having cross-sectional areas 
A as design variables (i.e., x° = v = A) are given by 
g*.(x) = Omax - o-(n l ) > 0. These constraints depend 
on the displacements x l = u\ and not on x° = A; so 
there is no z° that will satisfy the inequality in relation 
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by: 


(10) . Thus, a different decomposition is required. 

The formulation SDEP1 utilizes the second decom- 
position method described previously to enable satis- 
faction of the constraint qualification relations. Since 
functions for calculating the stress and buckling con- 
straints can be constructed by using both the sizing 
variables and the stress resultants of a structure, a 
change of variable is made to utilize stress resultants in 
the constraint functions. Because the stress and buck- 
ling constraint functions depend explicitly on the sizing 
variables with this change of variable, the satisfaction 
of the constraint qualification is much more likely, but 
satisfaction is not guaranteed as will be explained sub- 
sequently. The stress resultants within every element 
in the structure are computed within the structural- 
response subproblems, and are formally represented in 
subproblem i by the equation 

N* = N*(v, u*) (11) 

In formulation SDEP1, the structural- response sub- 
problem i consists of a solution of the unconstrained 
minimization given by statements (3) or (4) for u l fol- 
lowed by calculation of stress resultants N l by equation 

(11) . The design variables of structural- response sub- 
problem i are defined as x l = (u\ N’). The stress resul- 
tants N* are used instead of displacements in formulat- 
ing the stress and buckling constraints of the structural- 
sizing subproblem, which is also denoted subproblem 0. 
Thus, the structural-sizing subproblem in equilibrium 
programming formulation SDEP1 is: 

min WYv) 

V 

V - V min > 0 (12) 

g^C^N^N*)) > 0 

where i = The design variables of this 

structural-sizing subproblem are simply x° = v. 

The method of fully stressed design (ref. 5) can 
be derived from SDEP1 if: 1) only one-dimensional 
rod and two-dimensional membrane finite elements are 
used; 2) one sizing variable is associated with each fi- 
nite element having a stress constraint; 3) the stress 
constraints limit the maximum stress magnitude or von 
Mises stress; and 4) there are no buckling constraints. 
The structural- si zing subproblem can then be decom- 
posed into a set of independent elemental problems, one 
for each sizing variable and constrained element com- 
bination. The solution of these elemental problems is 
simple since the value of the sizing variable that makes 
a constraint active can be found analytically for each 
load case. For example, the elemental problem for cross- 
sectional area of one-dimensional rod element j is solved 


Aj — max(/lj |m i n , \P%jl\/<7j 7 m&x> ■ • • » I I A 7 ;, max ) 

(13) 

where the function max(...) chooses the maximum of 
its arguments, the structural-sizing design variable x® is 

defined to be the sizing variable Aj , and Pfj is the ax- 
ial force in the element j for load case i. The quantities 

P?f> are elements of the vectors N l which are also ele- 
* 

ments of x* and x°. A solution method that alternates 
between solving the sizing elemental problems, and the 
structural-response subproblems leads to fully stressed 
design. 

As stated previously, the change of variables that 
recasts the stress and buckling constraints in terms of 
sizing variables and stress resultants makes the satis- 
faction of the constraint qualification much more likely, 
but not guaranteed. A simple example makes this state- 
ment clear. Assume a rod is fixed between two rigid 
walls and its temperature is increased; the design prob- 
lem is to size the rod cross-sectional area A to minimize 
weight and satisfy a maximum stress constraint. The 
temperature change induces a strain which can be ex- 
pressed as an equivalent external load that is a func- 
tion of the stiffness. Thus, the equivalent external load 
and the rod stress resultant are both calculated in the 
structural- response subproblem as P z — kA AT where 
it is a constant of proportionality, and AT is the tem- 
perature rise. The stress, which is calculated in the 
structural-sizing subproblem, is given by a — P z /A , and 
also equals the quantity k AT. If this quantity is larger 
than the stress limit, it is intuitively obvious that no so- 
lution can be found, and in practice, the rod area using 
the SDEP1 formulation would increase without bound. 
The failure of SDEP1 to find a solution is a direct re- 
sult of the requirements for existence of an equilibrium 
point. One requirement is the boundedness of the feasi- 
ble region of the design space. This requirement can be 
satisfied by arbitrarily imposing a maximum area con- 
straint on the rod. However, at the feasible point where 
both the maximum area and the maximum stress con- 
straints are active, it is not possible to change the rod 
area to move into the interior of the feasible region, vio- 
lating the constraint qualification required for existence 
of an equilibrium point. Thus, the failure of SDEP1 
(and fully stressed design) to find a solution to this sim- 
ple problem is clarified by the EP existence theorem. 

Because SDEP1 is a form of fully stressed design, 
it shares the advantages and disadvantages of fully 
stressed design. The primary advantage is the simple 
nature of a structural-sizing subproblem that requires 
no derivatives and is easily decomposed into a set of in- 
dependent elemental problems. The disadvantages are: 
1) there is no mechanism to ensure satisfaction of the 


7 


necessary conditions represented by statement (6) so 
that the resulting equilibrium point may not be an op- 
timal point; and 2) constraints that have no explicit de- 
pendence on the sizing variables, such as displacement 
constraints, are not considered. An EP formulation is 
desired which satisfies the necessary conditions repre- 
sented by statement (6) at an equilibrium point, and 
can satisfy displacement constraints. The two formula- 
tions described subsequently overcome these disadvan- 
tages of SDEP1. 

Struct ura/- Design Formulation SDEP2 . As dis- 

cussed for SDEP1, if EP subproblems i = 1, . . . , M are 
identified as structural-response subproblems each rep- 
resented by statement (3), and subproblem 0 is identi- 
fied with a structural-sizing subproblem given by state- 
ment (5), the constraint qualification relations (10) can- 
not be satisfied for any displacement constraint because 
the displacements would be fixed parameters within 
the structural- si zing subproblem. In addition, the con- 
straint qualification relations may not be satisfied for 
certain stress and buckling constraints. Utilizing the 
third decomposition method described previously can 
overcome this difficulty. 

In EP formulation SDEP2, the displacements u\ 
which are parameters independent of the sizing vari- 
ables in structural-sizing constraint functions g(v,u l ), 
are replaced with a first-order Taylor series approxima- 
tion given by 


u 


A i 


(v,x’) = ^-(v- Vi) + 11* (14) 


In equation (14), the matrix du l /dv is a matrix of 
optimal sensitivity derivatives (ref. 8) of the displace- 
ments in subproblem i with respect to v, and v* is the 
value of v utilized in subproblem t when the sensitivities 
are calculated. All the design variables for structural- 
response subproblem i are utilized in equation (14) since 
x* = (u l ,duVdv, Vi). This approximation satisfies the 
following two properties. First, the approximation de- 
pends explicitly on v so that constraint qualification 
relations (10) necessary for equilibrium point existence 
can be satisfied. Second, it satisfies the conditions: 

* = u* and du Ai /dv = du'fdw at the equilibrium 
point x = x* where v = v* = v*. This second set of 
conditions ensures the optimality of the design given by 
the equilibrium point because the EP necessary condi- 
tions represented by statement (9) are then the same 
as the NLP necessary conditions represented by state- 
ment (6). 

Using the definition of equation (14), the structural- 


sizing subproblem is given by the following statement: 


min W(v) 

x°=v 

g*(v) > 0 

gufa^v.x’)) > 0 


(15) 


The structural-response subproblems for SDEP2 are 
the subproblems represented by unconstrained mini- 
mizations given by statement (3) (or the necessary con- 
ditions given by equation (4)), and the following equa- 
tions that determine the behavior variables du l /dv and 


v»: 


K< + ^ = [0| 


(16) 


Vi 


= V 


for i — l,...,Af. The first of the two equations in 
(16) can be viewed as an optimal sensitivity analysis 
of u* with respect to parameters v since u* are design 
variables of the NLP problem given by statement (3), 
and the second equation is a trivial identity which saves 
the values of the structural-sizing variables at which the 
sensitivity analysis was performed. 

SDEP2 is equivalent to a NLP approach to struc- 
tural design with a first-order Taylor series for an ap- 
proximate analysis when the equilibrium point is found 
by the following sequential steps: 1) solve the subprob- 
lems repr^ented by statement (3) (or necessary con- 
ditions represented by statement (4)) for u 1 for each 
i = l,...,Af; 2) solve equations (16) for du x /dv and 
3) use the quantities found in steps 1 and 2 in 
equation (14) to form approximate models for displace- 
ments u A *(v,x l ); 4) utilize the approximate models of 
step 3 in the structural-sizing subproblem represented 
by statement (15) to solve for v; and 5) repeat steps 
1 through 4 until the changes in the solutions from the 
previous iteration are smaller than a specified tolerance. 

Structural- Design Formulation SDEP3 . In the pre- 
ceding decompositions derived within the EP frame- 
work, all the design constraints are included in the 
structural-sizing subproblem. In the SDEP3 decom- 
position, the displacement constraints g{, are applied 
within structural- response subproblem t, so that the 
unconstrained structural-response subproblem, repre- 
sented by statement (3), becomes a constrained min- 
imization. This step corresponds to the first decom- 
position method described previously. However, the 
constrained structural-response subproblem no longer 
gives the response of the structure to the applied nodal 
forces F l E because “virtual” nodal forces are also ap- 
plied to satisfy the displacement constraints. These vir- 
tual nodal forces and the steps required to obtain the 
correct structural response are discussed subsequently. 
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The constrained minimization solved within struc- 
tural-response subproblem t is 


min Q(u < ) t K(v)u < - (F* £ ) T u < ) 

g'uK) > o 


(17) 


and the necessary conditions for this subproblem are 

K(v)u i -F , £ -(^) V = 0 

g’uCu 4 ) >0 (18) 

V > 0 


If the stress and buckling constraints are included 
within the structural-sizing subproblem, and they are 
calculated from stress resultants N’ of structural- 
response subproblem t, the second decomposition meth- 
od previously described is also utilized in SDEP3. As 
discussed for the SDEP1 formulation, this approach 
may not yield an optimum design for a redundant struc- 
ture. However, substituting approximate models, which 
depend explicitly on the sizing variables, for the stress 
resultants (as in reference 9, or equivalently, in the third 
method for decomposition) yields an optimal EP formu- 
lation. Thus, the stress resultants, and their sensitivity 
derivatives with respect to the sizing variables are re- 
quired, and are obtained from the equations 


(■V) T gu(u*) = 0 


It is apparent from the necessary conditions represented 
by statement (18) that quantities (<9g l u (u l )/<9u’) r A 1 
are similar to the applied nodal forces F^, and can 
be considered as “virtual” nodal forces. If the La- 
grange multipliers A 1 in these virtual nodal forces can 
be forced to zero by modification of the sizing variables, 
then the correct structural response is obtained. To ac- 
complish this task, the values of the Lagrange multi- 
pliers and their sensitivity derivatives with respect to 
the sizing variables must be supplied to the structural- 
sizing subproblem. Since statement (18) implys that 
the Lagrange multipliers and their sensitivity deriva- 
tives are zero for inactive constraints, and the sensi- 
tivity derivatives of the active displacement constraints 
are also zero, the sensitivity derivative equations to be 
solved in structural- response subproblem i are given by 
reference 8 as 


d 2 L’(u\ v, A*) du 1 
d(u l ) T <9u* dv 



^LV.y.A*) 

9(u 1 ) T, <9v 

dgu*<y) 

du* dv 


-[0] 


= [ 0 ] 


(19) 


In equations (19), the vector gu*(u l ) represents the 
active displacement constraints, the vector A 1 * are the 
Lagrange multipliers corresponding to these active dis- 
placement constraints, and the function L’(u’,v,A l ) is 
the Lagrangian function defined by 


£*(u*,v, A*) = l(u*) r K(v)u* - (F‘ £ ) r u < 

-(V)Vu( V) 

= 1(u*) t K(v)u* - (F* £ ) T u i 

- (K) T &\ L.(u’) 


( 20 ) 


r>T - N’(v> u‘) = 0 

tfl sp dN^v, u») du* dN^v.u*) ( 21 ) 

dv du’ dv dv 

In summary, structural-response subproblem t is 
represented by statements (17) or (18), and statements 
(19) and (21) with the addition of the trivial conditions 
Vi = v, as in method SDEP2, and the equations 


g‘ u = g' u ( ui ) 

d& = dg'uW du* (22) 

dv du* dv 

The design variables for this structural-response sub- 
problem are x’ — (u l , N\ A', g^, v*, du’/dv, dN’/dv, 
dXydv, dg,'Jd\). 

The formulation of the structural-sizing subprob- 
lem is now discussed. As mentioned previously, if the 
stress and buckling constraints are included within the 
structural-sizing subproblem and they are calculated 
from stress resultants N* of structural-response sub- 
problem t, then the substitution of approximate models 
which depend explicitly on the sizing variables for the 
stress resultants yields an optimal EP formulation. In 
SDEP3, a first-order Taylor series approximate model 
for stress resultant N’ is formed as 

*(v, x*) s (v - v<) + N* (23) 
dv 

where the matrix dN’/dv, defined previously, can be 
viewed as the total derivative of the stress resultants 
with respect to v. The quantities dN’/dv, Vj, and N’ 
are determined within structural-response subproblem 
i and are only parameters within the structural-sizing 
subproblem. 

Although the displacement constraints are satisfied 
within the structural- response subproblems, the La- 
grange multipliers corresponding to the active displace- 
ment constraints must be zero for the correct structural 
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response to be obtained. If a displacement constraint 
is not active within a structural- response subproblem, 
a provision within the structural-sizing subproblem is 
desired to guard against violation of the constraint. A 
constraint function that accomplishes this task within 
the structural-sizing subproblem is a first-order approx- 
imation in v to the expression g^ - A*. Component j 
of this approximation is given by: 

rf/(v,x<)={ 9 "> + ff <V ~ Vi) <9 “’ >0) (24) 

l ~ X j V ‘) (Sii; =°) 

Thus, if a displacement constraint is not active in 
structural- response subproblem i, the constraint func- 
tion in equation (24) is a fully linearized version of the 
corresponding form in the SDEP2 formulation. And, 
if a displacement constraint is active in structural- 
response subproblem t, then a structural-sizing con- 
straint with the functional form of equation (24) will 
coerce the virtual force in the structural- response sub- 
problem to zero. The constraint function in equation 
(24) satisfies all the requirements of the EP existence 
theorem described previously except for the continu- 
ity requirement. Continuity, and therefore equilibrium 
point existence, could be established if the function in 
statement (24) is replaced with one which differs from it 
only over the small transition interval 0 < ^ ■ < e. In 
this interval, the constant and gradient terms in equa- 
tion (24) are replaced with continuous transition func- 
tions that range between -A* and and between 

-dXj/dv and dg l u ^/dv , respectively. This continuous 
form of the constraint function satisfies the existence 
requirement, but has not been found necessary in the 
computational procedure. Therefore, the form of state- 
ment (24) is retained. 

With the previously mentioned definitions, the 
structural-sizing subproblem in SDEP3 becomes 


min W(v) 

x°=v 

gs(v) > o 

gffK^NKN^v.x'))) > 0 
gu ’(v.x 4 ) > 0 


(25) 


where i = 1, . . . , Af, the stress and buckling constraints 
depend on the stress resultants as in the subprob- 
lem represented by statement (12), and the vector of 
Misplacement” constraint functions g A *(v,x*) has ele- 
ments given by equation (24). A proof of the optimality 
of the equilibrium point obtained from SDEP3 is given 
in the appendix. 


Implementation of SDEP algorithms 

For practical structural-design problems, an equi- 
librium-programming-based solution algorithm should 
utilize existing finite element structural analysis soft- 
ware. Existing finite element programs can easily solve 
equations (4) for structural response, but are not well 
suited to solving a constrained NLP problem as repre- 
sented by statement (17) (or equivalently, the necessary 
conditions as represented by statement (18)). Because 
of this limitation, the class of displacement constraints 
considered in the present implementations of solution 
algorithms is restricted to simple displacement limits 
such as gu(u l ) = Umax — u l . In the following sec- 
tions, solution procedures are described for solving the 
SDEP2 and SDEP3 structural-response, and structural- 
sizing subproblems. These procedures are used with the 
Engineering Analysis Language (EAL) structural anal- 
ysis code (ref. 10) to solve some sample design prob- 
lems. The sequential solution approach described previ- 
ously, in which the structural-sizing subproblem and the 
structural-response subproblems are alternately solved, 
is the method chosen for the present study. 

Structural-Response Subproblems . The formula- 
tion of the structural-response subproblem in SDEP2 
uses standard methods for computing the structural re- 
sponse and the structural- response sensitivity deriva- 
tives (see ref. 11). Thus, existing finite element 
codes are well-suited to the structural- response sub- 
problem formulation of SDEP2. However, formulation 
SDEP3 specifies inequality displacement constraints 
in the structural-response subproblems represented by 
statement (17). The restriction of the displacement 
constraints to simple bounds on displacements in the 
present study allows for an easily implemented active 
constraint set approach to the structural- response sub- 
problems in SDEP3. In this approach, a structural- 
response subproblem is solved in an iterative manner 
so that a solution to the necessary conditions repre- 
sented by statement (18) is obtained from a sequence 
of solutions to equation (4) that differ by the boundary 
conditions assumed when factoring the stiffness matrix 
K. Specifically, when a displacement resulting from the 
solution to equation (4) violates either its maximum 
or minimum bound, the boundary conditions for the 
structural-response subproblem are modified to have 
that displacement specified to be the displacement limit 
when factoring K for the next solution to equation 
(4). The static reactions corresponding to these speci- 
fied displacements are computed by the finite element 
code. These static reactions equal the displacement- 
constraint Lagrange multipliers of statement (18) if the 
specified displacement is a lower bound, or the nega- 
tive of the displacement-constraint Lagrange multipli- 
ers if the specified displacement is an upper bound (see 
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the appendix). If any Lagrange multiplier determined 
from these static reactions is negative, the correspond- 
ing specified displacement should be “freed” in factor- 
ing the stiffness matrix for the next solution to equation 
(4). Likewise, if any displacement exceeds its allowable 
bounds, the appropriate displacement bound is utilized 
within the boundary conditions in the next iteration as 
previously described. The correct active set of displace- 
ment constraints is determined, and the correct solution 
to the necessary conditions represented by statement 
(18) are obtained when no calculated displacement vi- 
olates its bounds, and no Lagrange multiplier is nega- 
tive. An additional advantage of this active constraint 
set approach to solving the 5DEP3 structural- response 
subproblems is that after u l and A' are determined, 
the sensitivity derivatives du l /dv and dX x /dv are deter- 
mined by standard structural-response sensitivity anal- 
ysis methods, and it is not necessary to solve equa- 
tion (19). However, solving a constrained structural- 
response subproblem is typically more expensive than 
solving an unconstrained one because of the necessity 
of refactoring the stiffness matrix and resolving the sys- 
tem of equations for every change in the assumed set 
of active constraints. For the present study, the EAL 
sensitivity derivative runstreams of reference 12 were 
modified to generate sensitivity derivative information 
for both SDEP2 and SDEP3. 

StructuraJ-Sizing Subproblem. In the present study, 
the structural-sizing NLP subproblems represented by 
statements (15) for SDEP2 and (25) for SDEP3 are 
replaced with linear programming approximations to 
these subproblems having move limits as additional side 
constraints on the sizing variables. Move limits are 
easily incorporated within the EP theory given in ref- 
erence 6, so the details are omitted here. However, 
the method used for move limit control in the present 
study differs from reference 6 in that the move limit 
side constraints are modified before each solution of 
the structural-sizing subproblem, not after convergence 
to the equilibrium point corresponding to those move 
limits. The move limits chosen allow a 30% relative 
change in the design variables initially. The relative 
change allowed is reduced by 10% before each solu- 
tion to the structural-sizing subproblem, but it is not 
reduced below 6%. This strategy for controlling the 
move limits is fairly crude and has been found to be 
too restrictive for some of the example problems. In 
cases for which a linear programming approximation 
to a structural-sizing subproblems has no feasible so- 
lution, the objective function is modified to include a 
constraint violation penalty using the method described 
in reference 11. Thus, the subproblem will essentially 
minimize the worst constraint violation when no feasible 
solution exists. Also, because the stress constraints are 
linearized in this approach, there is, in theory, no dif- 


ference between the stress constraints calculated from 
either displacements or stress resultants. The subrou- 
tine DDLPRS of reference 13 was incorporated within 
EAL to solve the linear programming problems. 

For a large structural model, there can be thousands 
of elements that must be examined for local stress and 
buckling constraints, and must have the optimal struc- 
tural sizes determined. To reduce the number of design 
variables and constraints in the structural- sizing sub- 
problem, a large finite element model can be partitioned 
into regions. Within each region, all the elements are 
sized by a few design variables. In addition, the local 
constraints within each region can be “lumped” using 
the Kresselmeier-Steinhauser cumulative constraint de- 
scribed in reference 14. These reduction methods are 
utilized only for the civil transport wing problem de- 
scribed subsequently in the results section of the present 
paper. 

Results 

Two sample problems, a simple ten-bar truss exam- 
ple and a complex high-speed transport wing example, 
are evaluated using the SDEP2 and SDEP3 methods de- 
scribed above. SDEP2 is designated a “conventional” 
method herein, since it is equivalent to conventional 
NLP formulations using a Taylor series approximation 
for displacements for evaluating the constraints. In the 
present study, only a single load condition is consid- 
ered in the structural response. The results presented 
consider minimum gauge and displacement constraints 
only for the ten-bar truss example problem, and min- 
imum gauge, displacement, local buckling, and stress 
constraints for the high-speed transport wing example 
problem. 

Ten-Bar Truss 

The minimum weight ten-bar truss problem is illus- 
trated in figure 3. The vertical and horizontal mem- 
bers are each 360 in. long. The material properties 
assumed are those for aluminum with a Young’s mod- 
ulus of 10 7 psi, a Poisson’s ratio of 0.3, and a den- 
sity of 0.1 lb m /in 3 . Two 100,000 lb loads are applied, 
and the upper displacement limits ffj, 82 , <53, and £4 
for the displacement constraints are shown in the fig- 
ure. Results for two design cases are presented. In 
case 1, the limits for the displacement constraints are 

= 84 = 2.0 in. No 62 or <53 limits are assumed for 
this case. In case 2, the assumed displacement limits are 
6 \ = 2.0 in., 82 = 0.75 in., £3 = 0.5 in., and 84 = 1.8 in. 
No stress or local buckling constraints are considered 
in either case. The bar cross-sectional areas are the de- 
sign variables for the structural-sizing subproblem. The 
minimum cross-sectional area allowed is 0.1 in 2 . For 
both cases, results are obtained using two values for 
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Figure 3. Schematic diagram of geometry, loads, and dis- 
placement limits for ten-bar truss example problem. 


Table 1. Final bar cross-sectional areas and weight for ten- 
bar truss, case 1, initial area = 25 in 2 


Bar 

Number 

Ref. 11 
Area, 
in 2 

SDEP2 

Area, 

in 2 

SDEP3 

Area, 

in 2 

1 

30.52 

30.95 

31.05 

2 

0.10 

0.10 

- 0.10 

3 

23.20 

23.26 

23.23 

4 

15.22 

14.73 

14.77 

5 

0.10 

0.10 

0.10 

6 

0.55 

0,10 

0.10 

7 

7.46 

8.63 

8.63 

8 

21.04 

20.78 

20.69 

9 

21.53 

20.95 

21.02 

10 

0.10 

0.10 

0.10 

Weight, 

lb m 

5060.9 

5062.1 

5064.7 


the initial cross-sectional area of all members, 10.0 in 2 
and 25.0 in 2 . In the solution process, eighty steps of 
the sequential iterative solution method using the move 
limit strategy described previously are utilized for both 
the SDEP2 and SDEP3 formulations. 

Design case 1 results . The iteration histories of the 
truss weight resulting from formulations SDEP2 (i.e., 
the “conventional” formulation) and SDEP3 for the 
case having two displacement constraints are shown in 
figure 4 for the 25 in 2 initial cross-sectional area. In 
this figure and in subsequent figures, the value of the 
ordinate is the value at the beginning of the iteration 
indicated on the abscissa; therefore the initial design 
conditions are denoted as the conditions at iteration 1. 
The results for weight in figure 4 are very similiar for 


the two SDEP formulations. The iteration histories for 
the individual design variables also show little difference 
between the two formulations. The slow convergence of 
the weight is due to the restrictive move limit strategy 
utilized in the present study. The move limits control 
how fast the areas of bars 2, 5, 6, and 10 can be reduced, 
and do not allow these areas to achieve minimum gauge 
until iteration 63. Good agreement of final results 
from the two SDEP formulations with those reported in 
reference 11 is indicated by the final bar cross-sectional 
areas and truss weights given in Table 1. Because of 
small oscillations in the iteration history results, the 
SDEP3 results in this table are for the next-to-last 
iteration to present SDEP2 and SDEP3 results that are 
“in phase”. The agreement between the results from 
the two EP formulations is slightly better than their 
agreement with reference 11, but overall, the differences 
are minor. 

The iteration histories of the truss weights for this 
design case with the 10 in 2 initial cross-sectional area 
are shown in figure 5. The structural-sizing subprob- 
lems for both SDEP2 and SDEP3 have no feasible solu- 
tions for the initial two iterations (i.e., for results cor- 
responding to abscissa values 2 and 3 in figure 5). In 
these iterations, the structural-sizing subproblem mini- 
mizes the infeasibility, and both methods give identical 
results since they are at the side-constraint boundaries 
determined by the move limits. The results from the 
two SDEP methods differ for the next few iterations 
but again agree well at the later iterations. Again, the 
slow convergence to the optimum design is due to the 
restrictive move limit strategy. 

Design case 2 results . The iteration histories of the 
truss weights resulting from formulations SDEP2 and 
SDEP3 for the design case having four displacement 
constraints are shown in figure 6 for the 25 in 2 initial 
cross-sectional area. Again, there is little difference be- 
tween the weight histories of the two methods shown in 
figure 6, or in the individual design variable histories. 
However, much larger differences in the iteration histo- 
ries of weight occur for the 10 in 2 initial cross-sectional 
area. The results shown in figure 7, indicate that weight 
decreases for the initial iteration of SDEP3 while the 
weight increases for SDEP2. This initial decrease in 
weight for SDEP3 adversely affects the convergence of 
the method, and is discussed in detail subsequently. 
A closer examination of the results from structural- 
sizing subproblems for SDEP2 and SDEP3 reveals that 
the structural-sizing subproblem for SDEP2 has no fea- 
sible solution for the initial three iterations, but the 
structural-sizing subproblem for SDEP3 has no feasible 
solution for the initial eight iterations. For these iter- 
ations, the structural-sizing subproblem minimizes the 
infeasibility. For the 10 in 2 initial cross-sectional area, 
all the displacement constraints are initially active for 
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Figure 4. Convergence histories for weight of ten- bar truss 
example problem, design case 1 with initial area for 
sizing variables of 25 in 2 . 

the SDEP2 formulation, while only the displacements 
limited by 6 i, <52, and 63 are initially active in the 
SDEP3 formulation with corresponding virtual forces 
in the 6,000 - 85,000 lb range. For these infeasible iter- 
ations, to reduce the structural-sizing-subproblem ob- 
jective function of weight plus a penalty proportional 
to the largest virtual force (which corresponds to the 
82 displacement limit), the initial changes to the sizing 
variables reduce all the member areas except for bars 
7 and 8. Decreasing these member areas decreases the 
stiffness of the structure and reduces the virtual force 
required to satisfy the 62 displacement limit. No fea- 
sible solution is found in the next few iterations of the 
SDEP3 structural- sizing subproblem, but in these iter- 
ations more constraints contribute to the penalty and 
weight increases as more sizing variables begin to move 
toward the optimal values. The iteration histories for 
cross-sectional areas of bars 1, 4, and 6 are shown in 
figure 8 to illustrate typical sizing variable results. The 
significantly different sizing variable results in the initial 
iterations combined with the restrictive move limit con- 
trol strategy severely reduces the speed of convergence 
of the SDEP3 method in these results. 

High-Speed Civil Transport Wing 

The second sample problem considered in this study 
is the structural sizing of the wing for a proposed high- 
speed civil transport concept described in reference 15. 
The details defining this structural-design problem are 
too numerous to list so only a summary of its features 
is presented. The finite element model used for the 
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Figure 5. Convergence histories for weight of ten-bar truss 
example problem, design case 1 with initial area for 
sizing variables of 10 in 2 . 



Iteration number 

Figure 6. Convergence histories for weight of ten-bar truss 
example problem, design case 2 with initial area for 
sizing variables of 25 in 2 . 

structural -response subproblem is shown in figure 9. 
The upper wing cover panels are removed in this figure 
to illustrate the rib and spar web arrangement. The 
cover panels are titanium honeycomb sandwich panels, 
and the shear webs are titanium sine-wave webs. The 
model is relatively detailed with 1728 nodes (10,144 
degrees of freedom) and 2447 elements. A single load 
condition is analyzed in the structural response, a 2.5g 
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Figure 7. Convergence histories for weight of ten-bar truss 
example problem, design case 2 with initial area for 
sizing variables of 10 in 2 . 



Figure 8. Convergence histories for selected bar areas of 
ten-bar truss example problem, design case 2 with initial 
area for sizing variables of 10 in 2 * 

balanced, symmetric supersonic pull-up maneuver. 

There are 41 sizing variables considered in the 
structural-sizing subproblem. These design variables in- 
clude honeycomb sandwich faces heet thicknesses, hon- 
eycomb sandwich core heights, and sine- wave web 
gauges. The model is partitioned so that each siz- 
ing variable sizes multiple elements. The constraints 
considered are minimum and maximum gauges for the 



Figure 9. Finite element model of high-speed civil transport 
wing with upper cover panels removed. 


facesheets and webs, minimum and maximum thick- 
nesses for the honeycomb core, maximum element 
stresses, local buckling due to compression and shear, 
and a 12- foot limit on maximum deflection for the 2.5g 
maneuver. The constraint “lumping’ 1 method using the 
Kresselmeier-Steinhauser function described previously 
is utilized to reduce the large number of possible stress 
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and local buckling constraints from several thousand to 
only 106. The value of the Kresselmeier-Steinhauser 
parameter used in the present study is 50. 

Because of the computational expense of the de- 
tailed structural model, the two equilibrium program- 
ming formulations were executed for only 16 iterations. 
The relative move limit for the last iteration was thus 
slightly greater than the 6% minimum value. The his- 
tory of structural weight for the SDEP2 and SDEP3 
methods is shown in figure 10. The results from the two 
methods are nearly indistinguishable, and the weight 
reaches its limiting value of about 41,800 lb m by the 
eleventh iteration. The individual constraint having the 
largest effect on the structural weight is the 12- foot dis- 
placement limit. The convergence history of this con- 
straint for the SDEP2 and SDEP3 design methods is 
shown in figure 11. The scales for the displacement con- 
straint in SDEP2 and the corresponding virtual force 
constraint in SDEP3 are chosen to show equivalent vi- 
olation at the beginning of iteration 3 because at this 
iteration both methods have the same set of sizing vari- 
ables. The figure shows that the virtual force constraint 
in SDEP3 and the normalized displacement constraint 
in SDEP2 follow similar paths to reach feasible design 
space, however although both constraints are related to 
the wing displacement limit, they cannot be directly 
compared. The two methods yield differences in the 
design history for some of the wing sizing variables as 
illustrated in figure 12. Although the weight histories 
shown previously are nearly identical, noticeable differ- 
ences in the sizing variable histories for the two methods 
occur. Sizing variables 24 and 25 correspond to thick- 
nesses of shear webs with the normalized ordinate value 

1.0 corresponding to a thickness of 0.1 in. Sizing vari- 
ables 18 and 34 are typical, normalized honeycomb core 
and facesheet thicknesses, with normalizing factors of 

5.0 in. and 0.2 in., respectively. Although some differ- 
ences appear in the sizing variable histories of the two 
SDEP methods with the third iteration, the differences 
essentially disappear by the last iteration. 

Concluding Remarks 

The present study is an initial investigation in 
the use of equilibrium programming (EP) as a frame- 
work for design decomposition. The advantages of 
performing decomposition within the EP framework 
are the availability of existence theorems and optimal- 
ity criteria for the solution to the decomposed prob- 
lem. Three decompositions of structural- design prob- 
lems into structural-sizing and structural -response sub- 
problems are developed using this equilibrium program- 
ming framework. Two of the methods developed (i.e., 
SDEP1 and SDEP2) are shown to be equivalent to the 
existing conventional methods of fully stressed design, 
and of nonlinear programming using approximation 



Figure 10. Convergence histories for weight of high-speed 
civil transport wing example problem. 
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Figure 11. Convergence histories for displacement constraint 
of high-speed civil transport wing example problem. 


concepts for rapid analysis. In the third method for op- 
timal structural design presented (i.e., SDEP3), the dis- 
placement constraints are applied within the structural- 
response subproblems. However, the structural-sizing 
subproblem acquires new constraints which effectively 
control the magnitude of “virtual” forces that are ap- 
plied in the structural- response subproblems to sat- 
isfy the displacement constraints. The optimality of 
solutions generated by this third method is demon- 
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Figure 12. Convergence histories for selected sizing variables 
of high-speed civil transport wing example problem. 


strated analytically, and by numerical examples of a 
simple ten-bar truss design problem, and of a com- 
plex, high-speed civil transport wing design problem. 
In the example problems, the number of iterations to 
convergence for the third design method are approx- 
imately equal to the conventional nonlinear program- 
ming approach using first-order Taylor series approxi- 
mations for a rapid analysis. However, the total com- 
putational expense is slightly larger because of the ac- 
tive constraint set approach to solving the constrained 
structural- response subproblems. In addition, the con- 
vergence of the third method sometimes degrades when 
iterations occur which have no feasible solution to the 
structural-sizing subproblem. This degradation appears 
related to the tendency of the solution algorithm for 
the structural-sizing subproblem to minimize the mag- 
nitude of the “virtual” force constraint violation by in- 
creasing the flexibility of the structure. 


simplicity, only displacement constraints, and a single 
structural response corresponding to one load case are 
considered. Also, the displacement constraints are as- 
sumed to have the form 


gu = U max -U>0 (41) 

where it is also assumed that the nodal reference frames 
are oriented so that all deflection limits are maximums. 

The nonlinear programming formulation, such as 
given by statement (5), simplifies to: 

min W^v) 

v (42) 

gu = Umax - u(v) > 0 


where u(v) is the solution to the structural-response 
problem given by statement (3) with necessary condi- 
tions given by equation (4). In this nonlinear program- 
ming approach, the necessary conditions for statement 
( A2) are given by statement (6), two conditions of which 
simplify to 


dW(v) 

dv 


+ A ^=0 

ctv 

A* > 0 


(A3) 


where u* (and the corresponding A*) indicates those 
displacements which are at their constraint boundaries 
in the problem represented by statement (A2). The 
term du*/dv in the problem represented by statement 
(A3) is determined by taking the derivative of equation 
(4) with respect to v 


. du 9K(v)u r _. 

K(v) s; + -7sr L = |01 


(M) 


and forming the matrix du*/dv from terms of the 
matrix du/dv corresponding to the displacements at 
their constraint boundaries. 

In the SDEP3 equilibrium programming formulation 
of this problem, only those displacements at the con- 
straint boundary effectively contribute to the necessary 
conditions so the structural-sizing subproblem given by 
statement (25) simplifies to: 


Appendix 

The equivalence of the necessary conditions for the 
SDEP3 equilibrium programming method to a nonlin- 
ear programming approach that yields an optimal struc- 
tural design is demonstrated in this appendix. For 


dAj 

dv 


min WYv) 

x°=v v 

(v — vj) — Aj > 0 


(45) 


and the structural-response subproblem given by state- 
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ments (17) and (19) simplifies to: 


dition u = u l , yields the equation 


min 

x r 


in Q(«i 1 ) 7 ’K(v)u 1 -F£u 1 ) 


K ^ + 


[ 0 ] 
d \ I 

dv 


Umax — U* > 0 

V j _ v — 0 

dK{v)u l 


+ 


dv 


= [ 0 ] 


f?-n 


m 


K„ 

K 0 

( 

' du u " 
dv 



L 

' to] ' 


k d 

[ 

du» 

- dv - 


. [0] . 

) 

d\[ 

- dv - 


(49) 


where the matrix K has been partitioned as shown. The 
first set of equations in (A9) can be solved to give 



dUu -Lur-iif du * 
= -fr +K A K b 


(-410) 


where the formal definition of x 1 is x l = (u 1 , 

A*, v\, du l /dv, dA*/dv), u 1 is ordered into those dis- 
placements not at the constraint boundary u£, fol- 
lowed by those at the constraint boundary so that 
u 1 = (u£,uj), and the quantities \\ are the Lagrange 
multipliers for the displacements uj at the constraint 
boundaries. The two relations of interest in the neces- 
sary conditions for subproblem (A5) becomes 


dW(v) 

dv 


T dX\ „ 

+ " *r-° 

n > 0 


(•47) 


Substituting equations (A 10) in the second set of equa- 
tions in (A9) results in 

(k d - KlK- A 'K B ) £ = £ (/111) 

Finally, substituting equation (All) into the necessary 
conditions represented by statement (A7), the desired 
result of the necessary conditions represented by state- 
ment (A3) is obtained with 

Ar = /i r (K D -K^ 1 K 0 ) (A 12) 


where /i are the Lagrange multipliers for the inequality 
constraints in statement (A5). 

The desired result is to show that if x° = v and x 1 
is an equilibrium point for SDEP3 (i.e., solves (A5) and 
(A6) and thus their necessary conditions), then v and 
u = u 1 satisfy the necessary conditions for an optimal 
solution of the nonlinear programming formulation, in 
particular equation (4) and statement (A3). To accom- 
plish this result, express the two necessary conditions 
for the structural- response subproblem represented by 
statement (A6) as 


K(v)u 1 + 



-f e = o 

aJ > o 


(48) 


where the Lagrange multipliers corresponding to un- 
constrained displacements are explicitly shown as zero. 
Comparing necessary conditions represented by state- 
ments (4) and (A8), it is apparent that u cannot equal 
u 1 unless X\ = 0. This condition is assured by the 
second relation in statement (A8), and the inequality 
constraints in the structural-sizing subproblem repre- 
sented by statement (A5). Combining equation (A4), 
the last two equations in statement (A6), and the con- 


With the demonstration that the equilibrium point 
of SDEP3 satisfies equation (4) and the first relation in 
(A3), the final step required is proof that the inequal- 
ity in (A3) is also satisfied. To do so requires the as- 
sumption that the SDEP3 equilibrium point is isolated 
(i.e., there is no other equilibrium point within an c 
neighborhood of v). The proof will be by contradiction. 
First assume that at the equilibrium point of a SDEP3 
problem, one or more components of A* computed from 
equation (A 12) are negative. This assumption implies 
the existence of a more optimal, neighboring, feasible 
point. Specifically, within an e neighborhood of v and 
u, quantities v £ and u c exist and satisfy equation (4), 
the inequality constraint in statement (A2), and the 
relation W(v t ) < W(v). Since u c is the solution to 
an unconstrained potential energy minimization (i.e., 
equation (4)) with the stiffness matrix K(v e ), and it 
also satisfies the displacement constraints in statement 
(A2), then u 1 = u* is the solution to the constrained 
potential energy minimization given in statement (A6) 
for the same stiffness matrix. Thus, the necessary con- 
dition equations in both statements (4) and (A8) must 
be satisfied, and Aj = 0 for both equations to hold. 
Finally, the choice x° = vj = v e satisfies the inequal- 
ity constraint in statement (A5) with the weight W in 
statement (A5) smaller for x° = v t than for x° = v. 
Since the e chosen is arbitrary, either the original point 
is not an equilibrium point or the equilibrium point 
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is not isolated, contradicting the original assumptions. 
Thus, the inequality in statement (A3) must hold, and 
the equilibrium point is also an optimal point. 
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